Tunable capacitive inter-dot coupling in a bilayer graphene double quantum dot 
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We report on a double quantum dot which is formed in a width-modulated etched bilayer graphene 
nanoribbon. A number of lateral graphene gates enable us to tune the quantum dot energy levels 
and the tunneling barriers of the device over a wide energy range. Charge stability diagrams and 
in particular individual triple point pairs allow to study the tunable capacitive inter-dot coupling 
energy as well as the spectrum of the electronic excited states on a number of individual triple 
points. We extract a mutual capacitive inter-dot coupling in the range of 2 - 6 meV and an inter-dot 
tunnel coupling on the order of 1.5 /xeV. 



I. INTRODUCTION 



Graphene is a promising material for future quan- 
tum devices and quantum information technology 
In particular, graphene quantum dot devices are inter- 
esting candidates for spin qubits with long coherence 
times. The weak spin-orbit coupling Q and weak hy- 
perfine interaction is expected to be of advantage com- 
pared to state-of-the-art GaAs-devices. However, the 
gapless band structure of graphene makes it hard to elec- 
trostatically confine charge carriers. This can be over- 
come by nanostructuring graphene, where it has been 
shown that an etching-based "paper cutting" technique 
leads to graphene nanodevices where transport is domi- 
nated by a disorder-induced energy gap. Consequently, 
graphene nanoribbons single-electron transistors 

(SETs) [iMl], graphene quantum dots (QD) 
even with charge detectors [3: [3 H^l and graphene dou- 
ble quantum dots [2T1425j have been demonstrated suc- 
cessfully. These devices allowed the experimental obser- 
vation of excited states |22|. [26| , spin states [l^ and the 
very few carrier regime [2^. Almost all of these stud- 
ies were based on single-layer graphene and showed de- 
vice limitations related to the presence of edge rough- 
ness, disorder and unintended vibrational degrees of free- 
dom [29| . Bilayer graphene is an interesting candidate 
to overcome some of these shortcomings. In particular, 
bilayer graphene allows to open a band gap by an out- 
of-plane electric field [10, [IJ , which may enable a soft 
confinement potential and may reduce the influence of lo- 
calized edge states. Furthermore, it is expected that rip- 
ples and substrate-induced disorder are reduced, which 
may also increase the mechanical stability suppressing 
unwanted vibrational modes. 

Here we present a bilayer graphene double quantum dot 
with a number of lateral gates. The gates allow to tune 
transport over a wide energy range and they enable to op- 
erate the device in different configurations. We focus on 
the double quantum dot (DQD) configuration and show 
characteristic honeycomb-like charge stability diagrams. 
More specifically, we discuss the electrostatic tunability 
of the mutual capacitive inter-dot coupling energy. Fi- 
nally, we also extract an inter-dot tunnel coupling en- 
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FIG. 1. (a) Scanning force microscope image of the mea- 
sured device. The dashed lines highlight the contours of 
the graphene (light gray), (b) Raman spectrum recorded on 
the used bilayer graphene flake. The inset shows an opti- 
cal microscope image of the contacted bilayer graphene de- 
vice (the arrow points to the quantum device and the white 
scale bar is 20 /im) . (c) Back gate characteristics recorded 
at Vb = — 20mV. (d) Source-drain current as a function of 
back gate voltage (Vbg) and voltage applied to the side gates 
(Vig=Kg=14g) at Vb = — 20mV (after charge rearrangement). 



ergy. 



II. DEVICE FABRICATION 



The device is fabricated by mechanical exfoliation of 
natural graphite by adhesive tapes [H, [11]. The sub- 
strate material consists of highly doped silicon (Si"'""'") 
bulk material covered with around 300 nm of silicon ox- 
ide (Si02) in order to identify bilayer graphene flakes due 
to absorption and interference. Before deposition, refer- 
ence alignment markers are patterned on the substrate 
to relocate the flakes for further processing and charac- 
terization. In particular Raman spectroscopy measure- 
ments are used to identify bilayer graphene flakes [3^[35j . 
The graphene flakes have to be structured to submicron 
dimensions in order to fulfill the device design require- 
ments. We use a technique based on resist spin coating, 
electron beam lithography (EBL), development and sub- 
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sequent etching of the unprotected graphene. We use 
an EBL resist [polymethylmethacrylate (PMMA)] with 
a thickness of 80 nm and short etching time to define 
the nanoscale structures. We follow earlier work (l3j . 
where it has been shown that short (5 s) mainly physi- 
cal reactive ion etching (RIE) based on argon and oxy- 
gen (80:20) provides good results without influencing the 
overall quality of the graphene flake [s^. After etch- 
ing and removing the residual EBL resist, the graphene 
nanostructures are contacted by an additional EBL step, 
followed by metalization and lift-off. Here we evaporated 
5 nm chrome (Cr) and 50 nm gold (Au) for contacting 
the bilayer graphene double quantum dot device. 
A scanning force microscope image of the investigated 
device is shown in Figure [ija), where the graphene struc- 
tures (bright areas) resting on Si02 (dark areas) are 
highlighted by the white dashed lines. The bilayer dou- 
ble quantum dot consists of three 30 nm wide and 50 to 
100 nm long constrictions, two are separating the source 
and drain leads from the double quantum dot and one is 
separating the left (LQD) from the right (RQD) quantum 
dot. Both bilayer graphene quantum dots have a diame- 
ter of approximately 50 nm. In this study we use the left 
(Ig) and the right (rg) gate to change the number of car- 
riers in the left and the right QD, while the central gate 
(eg) is used to tune the inter-dot coupling. Additionally, 
the back gate (bg) is used to adjust the overall Fermi 
level. The Raman spectrum shown in Figure 1(b) shows 
that the investigated device is indeed based on a bilayer 
graphene flake (see also Fig. 1(b) in Ref. [1^). The 
FWHM of the 2D-peak is 57.9 cm"! with a 18.8 cm"! 
spacing of the two inner Lorentzian peaks which in com- 
bination with the shoulder on the left side (see 2D peak 
in Figure 1(b)) provides a good finger print for the bi- 
layer nature [SJ,!!!]. The inset in Figure 1(b) shows an 
optical image of the contacted device (see arrow) . 



III. MEASUREMENTS 

The measurements have been performed in a dilu- 
tion refrigerator at a base temperature of 10 mK. We 
have measured the two-terminal conductance through 
the bilayer graphene double quantum dot by applying 
a symmetric dc bias voltage while measuring the current 
through the device with a resolution better than 50 fA. 



A. Device Characterization 

To characterize the double quantum dot device on a 
large energy scale the back gate voltage is tuned over a 
range of GOV. Figure 1(c) shows the source-drain current 
as a function of the back gate voltage. We observe a so- 
called transport gap [7| above Vbg > 50 V (Vb = — 20mV) 
and hole dominated transport for lower back gate volt- 
ages (see inset in Fig. 1(c)). The large-scale current fiuc- 
tuations in the transport gap region of the DQD device 
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FIG. 2. (a) Source-dram current as a function of the left gate 
(Vig) and right gate (Vrg) voltage at an applied bias of (Vb ~ 
8mV) (b) Illustration of a triple point pair. The quantities 
are used to deduce the coupling energy, (c) Capacitive inter- 
dot coupling energy as a function of the central gate voltage 
Vcg for left (open circle) and right (full circle) gate values. 
The insets show three corresponding triple point pairs. 



are due to local resonances in the bilayer graphene con- 
strictions acting as tunneling barriers. In Figure 1(d) we 
show a measurement of the current as a function of the 
back gate voltage and the left, right and central gate volt- 
ages {V\g = Vrg = Vcg) at constant bias (Vb = — 20mV). 
One characteristic slope is observed, which separates the 
area of suppressed current (upper right corner) from the 
area of elevated hole transport (lower left corner). The 
three tunneling barriers can be tuned almost equally by 
either the back gate or the sum of the left, right and cen- 
tral gates (a = AI4g / AVig^rg^cg = 1-5), which is in a 
good agreement with the device design and with earlier 
studies [111. 

By fixing the back gate and the central gate voltage we 
can record charge stability diagrams by varying the left 
and right gate voltages independently. In Figure 2(a) 
we show a finite bias (Vb = 8mV) measurement, where 
the honeycomb pattern is observed, which is character- 
istic for the charge stability diagram of a DQD [33|. It 
is crucial that electrical transport through the DQD is 
only possible in the case where the energy levels in both 
quantum dots are aligned within the source-drain bias 
window. From this measurement we can extract the lever 
arms and the single quantum dot addition energies, such 
as eI;^^ = 23meV and E^^f = 13mcV for the left and 
right quantum dot, respectively. 
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FIG. 3. (a) Mutual capacitive inter-dot coupling energy plotted as a function of central gate voltage. The corresponding 
bias and electron configuration (M, N) in the (left, right) quantum dot are given in the upper left corner. The inset shows 
the position of the triple points in the charge stability diagrams revealing (M, N). The corresponding color scale bar is shown 
in panel (c) and has to be multiplied by a factor of -0.5. (b,c) Current through the double dot as a function of the side gate 
voltages for different central gate voltages. The very same triple point is measured for Vb=6mV (b) and Vb ~ — 6mV (c), for 
increasing Vcg, ranging from 7.40 V to 7.58 V. The capacitive inter-dot coupling energies for (b) and (c) are represented by the 
black and gray line in panel (a). The black cross in panel (b) highlights the adjustment of the triple point position due to 
capacitive cross-talk. The dotted line in panel (c) highlights the position of a charge rearrangement. For each row the color 
scale bar is located on the very right. In the first three triple point pairs in panel (c) the current has to be multiplied by a 
factor 2. 



B. Capacitive Inter-dot Coupling 

In this section we discuss the effect of the central gate 
voltage on the capacitive coupling between the left and 
the right quantum dot (i.e. the capacitive inter-dot cou- 
pling). In Figure [2Ib) we show a schematic illustration 
of a triple point pair highlighting the voltage differences, 
which allow to determine the mutual capacitive coupling 
energies, i?™ and E^^. These coupling energies can be di- 



rectly extracted from the lever arms and the illustrated 
right gate and left gate voltage differences, respectively. 
For example, we obtain E"™ from 

where Q!rg,R = Vh/6Vi-g is the right gate lever arm. A sim- 
ilar expression exists also for left gate voltage differences 
exclusively. In Figure [2^c) the coupling energies ^-i^ig, 
obtained by the right gate voltage values (full circles) and 
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by the left gate voltage values (open circle), are plotted 
for eight different central gate voltages for one and the 
same triple point pair. The three insets show examples 
of the investigated triple point pair for a weak, medium 
and stronger inter-dot coupling. It is expected that E™ 
and have the same values, but due to limited res- 
olution and partial charge rearrangements the edges of 
the triangles are not very clear, leading to a scattering of 
both results. 

The average values i?™ = (i?™ -I- E^)/2 as function 
of the central gate voltage are shown in Figure IHJa) (see 
stars ir). The mutual capacitive coupling energy is in- 
creasing non-lincarly as function of the central gate volt- 
age. Next we study the tunability of the capacitive cou- 
pling energy for near by triple point pairs. Since the 
absolute number of electrons in the two quantum dots 
are not known, we focus on relative changes and assume 
as a starting point M electrons in the right and N elec- 
trons in the left dot. For a near by triple point pair 
with one additional electron in the right dot (M-|-1,N) 
also a positive slope for E™ versus V^g is observed (see 
diamonds ♦ in Fig. [31[a)). However, for the triple point 
pair with (M-|-5,N-1) electrons the capacitive inter-dot 
coupling energy is decreasing as function of the central 
gate voltage (see full squares ■ in Fig. |3Ua)). 
The corresponding triple point pairs from which the cou- 
pling energies have been extracted are shown in Fig.[3Jb). 
With increasing central gate voltage the spacing between 
the triangle tips is decreasing, which is directly related 
to the decreasing capacitive coupling energy. A closer 
inspection of the excited states shows not only a cou- 
pling energy dependent reduction of the current within 
the triangles, but also a change in the coupling of the 
excited state with an energy of w 2 meV (sec line run- 
ning parallel to the base-line; highlighted by the black 
arrow). In the first row of Figure [3jb) the excited state 
resonance is almost completely suppressed and becomes 
more clear when increasing the central gate voltage (see 
second row) . Due to capacitive cross-talk the overall po- 
sition of the triple point pair changes in the left and right 
gate voltage plane as function of the central gate voltage. 
This is highlighted by the distance to the fixed black cross 
(-I-), which is fixed to the very same position in the gate 
voltage plane. 

For a negative bias (Vb = —6 mV) the same dependence is 
observed (see open squares □ in Fig.[3][a) and Fig. [31Jc)). 
Although the absolute value of the bias is the same, the 
trace with the open squares has an offset, compared to 
trace with the full squares. Please note that the corre- 
sponding triple point pairs (see Fig. [3Jc)) nicely show 
the effect of the central gate on a charge rearrangement 
(see vertical dotted line). In the first two panels the 
rearrangement is located at the left edge of the triple 
point. With increasing central gate voltage the triple 
point pair is moving to lower right and left gate volt- 
age values. For the rearrangement, which is located very 
close to the right barrier (vertical dotted line) , the cross- 
talk is significantly lower and therefore it moves through 
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FIG. 4. (a) Line cuts along the detuning energy axis, as 
shown in panel (b). The bold lines are Lorentzians fitted to 
the ground state, (b) Current as a function of the side gate 
voltage Kg and Vig for an applied bias of Vt, = 6 mV and Vb 
= -6mV, respectively. 



the triangles and lies outside for central gate voltages 
above Kg > 7.58 V. 

Overall, we observe coupling energies ranging from about 
2 meV to 6 meV. This is by a factor of 2 larger, compared 
to values from an earlier study on a single-layer double 
quantum dot plj . This difference might be explained by 
the considerably smaller size of the investigated quan- 
tum dots in our work. More importantly we also observe 
a rather non-monotonic behavior of the capacitive cou- 
pling strength as function of Vcg and (M,N). This is most 
likely due to a disorder dominated electrostatic landscape 
in the narrow constriction connecting the two quantum 
dots, which is in agreement with recent experiments on 
graphcnc tunneling barriers jToj . 



C. Inter-dot Tunnel Coupling 

Transport through the quantum dot states can be ana- 
lyzed more quantitatively using the result from Stoof and 
Nazarov for resonant tunneling [38| . In the limit of weak 
intcr-dot tunnel coupling ^ Fin^out, where tm is the 
inter-dot tunnel coupling, and Fin. out are incoming and 
outgoing tunnel rates, the current / follows a Lorentzian 
line shape as a function of the detuning energy e, 



1 + {2e/hTa 



(3.2) 



where h is the Plank constant. In Figure 4(a) we plot 
two line-cuts along the detuning axis for a positive and 
a negative bias (Vb — i 6mV), and Lorentzian line fits 
to the ground state lines. Please note that in partic- 
ular in the lower panel of Fig. 4(a) the excited states 
at an energy of 1.8 meV and 3.6 me V are nicely visible 
(see arrows). The corresponding triple points with con- 
stant central gate voltages of V^g = 7.465 V (upper panel) 
and V^g ~ 2Y (lower panel) are shown in Figure l^^b). 
Following Ref. [23 the fitting is done only for the data 
points outside of the bias triangle in order to minimize 
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the contributions from the inelastic transport. We ex- 
tract from the fittings a tunnel rate from the right dot 
to drain /irR_j.d « 210 /leV and an inter-dot tunnel rate 
htm ~ l.bfieV. From negative bias we obtain a tunnel 
rate from the left dot to source /iFl-s-s ^ 420 fieV and an 
inter-dot tunnel rate htm ~ 1.7//eV. 
These values are in reasonable agreement with the results 
reported by Liu et al. [2^. The out of the dot tunnel 
rates have the same order of magnitude. The inter-dot 
tunnel rate, however, is in our case by a factor of six 
smaller (<m,Liu = 10/j,eV), which might be explained by 
the different device design. In their case the double dot 
is formed by two top gates on a 20 nm wide graphene 
nanoribbon. The tunneling barriers, induced and con- 
trolled by the 40 nm wide central gate, seem to be more 
transparent, compared to our side gated 30 nm wide con- 
striction. Please note that the extracted inter-dot tunnel 
rate is also consistent with the upper bound of t < 20 ^eV 
estimated by Molitor et al. [22!| for a single-layer graphene 
double quantum dot. 

IV. CONCLUSION 

We have fabricated a bilayer graphene double quan- 
tum dot device, which can be tuned by a number of 
all-graphene lateral gates. The device has been char- 
acterized on a large energy range and its functionality 



is demonstrated by charge stability diagrams and well- 
resolved triple points. 

We have studied the mutual capacitive inter-dot coupling 
energy E"^ and estimated the tunnel coupling energy tm- 
By varying the central gate voltage we are able to tune 
E"^ systematically from 2 to 6meV which is in good 
agreement with former results for single-layer graphene 
double quantum dots (23 . [23| . taking into account the 
smaller size of the presented quantum dots. Further- 
more, the detailed investigation of individual triple point 
pairs has allowed to determine electronic excited state 
energies and an inter-dot tunnel rate of around 1.5 ^eV 
which is roughly one order of magnitude lower, compared 
to recent results reported for single-layer graphene double 
dots [12, . Most likely this difference might be due to 
the different device designs and architectures. However, 
further investigations and more statistics are needed to 
elucidate this issue. 

In summary these results give insights into tunable bi- 
layer graphene double quantum dot devices and the abil- 
ity to tune the capacitive inter-dot coupling energy which 
is promising for the implementation of spin qubits in fu- 
ture graphene devices. 
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cussions with J. Glittinger and S. Droscher and support 
by the JARA Seed Fund and the DFG (SPP-1459 and 
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